Proc. Int. Conf. Heavy Electrons (ICHE2010) 
J. Phys. Soc. Jpn. 80 (2011) SA134 
© 201 1 The Physical Society of Japan 



Electron Mass Enhancement due to Anharmonic Local Phonons 

Kunihiro Oshiba and Takashi Hotta 

Department of Physics, Tokyo Metropolitan University, Hachioji, Tokyo 192-0397, Japan 



o 

(N 

Oh- 
< 

(N 



u 

• 

a 

o 

(N 
> 

(N 

m 

(N 
O 



In order to understand how electron effective mass is enhanced by anharmonic local oscillation of an 
atom in a cage composed of other atoms, i.e., rattling, we analyze anharmonic Holstein model by using a 
Green's function method. Due to the evaluation of an electron mass enhancement factor Z, we find that Z 
becomes maximum when zero-point energy is comparable with potential height at which the amplitude of 
oscillation is rapidly enlarged. Cooperation of such quantum and rattling effects is considered to be a key 
issue to explain the electron mass enhancement in electron-rattling systems. 
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In Sm-based filled skutterudite SmOs4Sbi2/^ it has been 
observed that electronic specific heat coefficient 7e is signifi- 
cantly enhanced, suggesting the emergence of heavy-electron 
state. An intriguing issue is that 7e is almost unchanged even 
when we apply a magnetic field up to 30 Tesla. Namely, the 
heavy-electron state is magnetically robust. On the basis of 
the standard Kondo effect originating from spin degree of 
freedom, it is difficult to explain that heavy effective mass 
is independent of the applied magnetic field. Thus, the non- 
magnetic Kondo effect due to electron-phonon interaction has 
been pointed out,^^ since local anharmonic oscillation of rare- 
earth atom, i.e., rattling, in the pnictogen cage has played a 
crucial role for electronic properties of filled skutterudites. 

On such a background, Kondo effect with phonon origin 
and its related phenomena have been vigorously investigated. 
As an extension of the two-level Kondo problem, the four- and 
six-level Kondo systems have been analyzed to understand 
the effect of rattling in filled skutterudites.^' The origin of 
heavy mass has been discussed in periodic Anderson-Holstein 
model^^ and in Holstein model. Due to the evaluation of 7e 
under the magnetic field in Anderson-Holstein model, it has 
been found that large 7e becomes magnetically robust, when 
the bottom of the potential is relatively wide and flat. From 
these efforts, it has been gradually recognized that magnet- 
ically robust heavy-fermion phenomenon actually occurs in 
electron-rattling systems, but it is further necessary to clar- 
ify properties of the heavy-electron state in electron-rattling 
systems from various aspects. In particular, the analysis in the 
periodic system is required. 

In this paper, we analyze anharmonic Holstein model by 
using the standard Green's function method, in order to un- 
derstand the mechanism of electron mass enhancement due 
to rattling. For the purpose, we evaluate an electron mass en- 
hancement factor Z by focusing on the relation between zero- 
point energy Eq and a potential height Vq. Note that Vq is 
the potential energy at which the amplitude of oscillation is 
rapidly increased. For Vq > Eq, Z exhibits a peak aiT ^ Vq, 
where T is a temperature, indicating that the amplitude of os- 
cillation is relevant to the enhancement of Z. For Vq < Eq, 
we do not find the peak in Z, but the shoulder structure is ob- 
served at T ~ £^0- When we change the states between the 
cases of Vq > Eq and Vq < Eq, it is found that Z takes the 



maximum value for Vq ^ Eq. Thus, we envisage a scenario 
that Z is significantly enhanced when both quantum and rat- 
tling effects work cooperatively. 

Let us introduce anharmonic Holstein model, given by 



(1) 



where Cka indicates an annihilation operator of electron with 
momentum k and spin a, Sk denotes electron energy, g is an 
electron-phonon coupling constant, i denotes a site, Qi is nor- 
mal coordinate of oscillation of an atom at site i, and Cia- is 
an annihilation operator of electron at site i. Throughout this 
paper, we use such units as h=kB=l- 

The third term in eq. (1) indicates the oscillation of atom, 
given by i^ph=Xli [^i^/^ + ^iQi)]^ where Pi is the canonical 
momentum and V{Qi) is an anharmonic potential for atom. 
Note that the reduced mass of the oscillation is set as unity. 
The potential is expressed by V{Qi)=uj'^Qj/2-\-k4Qj-\- k^Q^, 
where uj is the energy of oscillator and k^ and k^ denote the 
coefficients for forth- and sixth-order anharmonic terms, re- 
spectively. Here we introduce the phonon operator defined 
through Qi=(ai + a])/ V2uj. Then, we obtain 

i/ph =c^ ^[alai + l/2 + /3(ai + al)^+7(a, + al)% (2) 

i 

where /3 and 7 are non-dimensional coefficients of anhar- 
monic terms, given by /3=k4 / (Aoo^) and ^=ke / {Suj"^) . By us- 
ing such non-dimensional parameters, we express the poten- 
tial as V{qi)=auj{qj + Wa^qj + 64a^7gf ), where a is the 
non-dimensional electron-phonon coupling constant given by 
OL-g^ I and qi-Qiji. Here i denotes the typical length 
for oscillation, given i-gjuJ^. 

First let us briefly discuss how the potential shape is 
changed by anharmonic parameters. Note that in this paper, 
we consider only the case of 7 > 0, since the oscillation 
should be confined in a finite space. When we change the 
value of /3, the potential shapes are classified into three:^^ 
(i) On-center type for — \/37/4 < /3 < 0, (ii) rattling type 
for — < < — Y^37/4, and (iii) off-center type for 
/3 < — y^. As shown in Fig. 1, the on-center type potential 
has a minimum at g=0, while in the off-center type, we ob- 
serve three minima. The rattling-type potential has relatively 
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Fig. 1. (Color online) Change of the potential shape at 7=10~^. Dotted 
curve denotes harmonic potential shape for comparison. The curve of 
/3=— 0.001 denotes the on-center type, while that of /3=— 0.0035 indicates 
the off-center one. The curves of /5=— ^^87/4 and are categorized 

into the rattling type. In this paper, we focus on the case of /3=— ^^87/4. 

wide region in the bottom. Since we are interested in the effect 
of oscillation in the rattling-type potential, here we consider 
the case of y37/4, depicted by thick curve in Fig. 1. In 
such a case, the potential has two saddle points and we de- 
fine Vb as the potential height at the saddle point. When we 
increase the energy of oscillator, the amplitude of oscillation 
is considered to be rapidly enlarged at the energy of . 

Now we calculate the electron self-energy S in the 
secound-order perturbation in terms of the coupling g within 
an adiabatic approximation. Namely, we consider a situa- 
tion in which electron bandwidth W is much larger than the 
phonon energy cj. Hereafter W is taken as the energy unit 
and we set W-X. First we evaluate the bare phonon Green's 
function D. For the purpose, we diagonalize i^ph to obtain 
eigenenergy En and eigenstates |n). Then, D is obtained in 
the Lehmann representation as 

^(")=oE '" \{m\{ai^a\)\n)\\ (3) 

i I — Z + hjn — -L^m 

where 2; is a complex number with the dimension of energy 
and ft is 3. partition function. Note that the site dependence 
does not appear in D, since we consider Einstein-type local 
phonons. Then, we evaluate the electron self-energy as 

T^{iuJn) = -auj'^T^^D{iuJn - iuJn^)G{k\iuJn^), (4) 

n' k' 

where G denotes the electron Green's function given by 
G{k^iujn)='X/ {iujn — Sk) and ujn is fermion Matsubara fre- 
quency defined by (jJn={2n + l)7rT with an integer n. By as- 
suming the electron density of states with rectangular shape of 
the width W, we perform the summation in terms of k\ Note 
that the half-filling case is considered in this paper. Then, we 
obtain E as 




T (unit of W) 

Fig. 2. (Color online) (a) Mass enhancement factor Z vs. temperature T 
for 7=0.0005 with Eq < Vq and for 7=0.0 07 wi th Eq > Vq. We note 
that we always keep the relation of /3=— ^87/4. Note also that Ti is 
a peak temperature, while T2 is a temperature at which Z begins to be 
constant with decreasing T. (b) Non-dimensional Debye- Waller factor / 
vs. T and (c) df /dT vs. T for the same parameters as (a). 

where A denotes the Eliashberg electron-phonon coupling 
constant, defined by X=2auj/W. In actual calculations for 
E, we exploit a Fast-Fourier-Transformation technique. Then, 
we define the mass enhancement factor Z as 

Z = 1- Y.{mT)/{mT). (6) 

Here we use the value on the imaginary axis at the lowest tem- 
perature as an approximation for Z, although Z should be de- 
fined by the energy differentiation of E on the real axis. Since 
E is found to exhibit a peak at the energy in the order of cj, we 
replace the energy differentiation with the energy difference 
of the width T iovT < uo. Thus, we adopt the present defini- 
tion for Z to discuss properties of Z in the region of T < a;. In 
fact, characteristic structure of Z is found for T < cj, which 
is consistent with the present condition. 

Now let us discuss the calculated results. Here we set 
A=1.0 and a;=0.1. For the diagonalization of i^ph, we use 250 
phonon basis. In the calculation of E, we take 4096 Matsub- 
ara frequencies. First we discuss temperature dependence of 
Z by focusing on the relation between Vb and zero-point en- 
ergy £^0- In Fig- 2(a), we depict two curves for 7=0.0005 and 
0.007. Note that the condition of y^37/4 is always kept. 
For 7=0.0005, we find that £:o=0.00432 and ^0=0.0108, sug- 
gesting the case of Vb > Eq. On the other hand, for 7=0.007, 
we obtain £;o=0.00359 and l/o=0.00288, leading to the situa- 
tion with Vq < Eq. 

As observed in Fig. 2(a), for the case of 7=0.0005, Z has a 
broad peak at a temperature which is defined as Ti in this pa- 
per. From the detailed calculations, it is found that Ti is scaled 
by Vq. When we further decrease the temperature, we observe 
the structure of a shoulder. We define T2 as the temperature at 
which Z becomes constant. For the case of 7=0.007, on the 
other hand, we find no peak and only T2 can be defined. Af- 
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ter performing the calculations for several parameter sets , we 
have found that T2 is scaled by Eq. It is also found that the 
first-order excitation energy is always in the order of £^0 in 
the present parameters, suggesting that the temparature scale 
is considered to be £^0 at low temparatures. 

From these results, first we claim that the enhancement of 
Z is closely related to the potential height Vq. By taking into 
account the definition of Vq, we point out that Z becomes 
large when the amplitude of oscillation is rapidly increased. 
This should be called the rattling effect. The result is consis- 
tent with the previous conclusion from the numerical calcu- 
lation of 7e in Holstein- Anderson model. Second we claim 
that the rattling effect is not significant for T < Eq, since 
zero-point oscillation is dominant in such a region and the 
rattling effect is masked. 

In order to examine the relation between Z and the am- 
plitude of oscillation, we evaluate the Debye- Waller factor 
which represents the intensity of thermal motion of the caged 
atom. In the non-dimensional form, it is given by 

^ \G\^P 3£2 ' ^'^ 

where Fdw denotes the Debye- Waller factor, G denotes the 
reciprocal lattice vector, and (• • • ) indicates the operation to 
take thermal average. Note again that the site dependence is 
ignored, since we consider Einstein-type local phonons. 

In Fig. 2(b), we show calculated results for /. In order to 
clarify the variation of / in terms of T, we also plot df /dT 
vs. T in Fig. 2(c). From these results, we claim that / rapidly 
increases at T Ti for the case of Vq > Eq. If we further 
decrease T, we arrive at a constant value at T2 ^ Eq. For the 
case of Vq < Eq, however, the value of df /dT is not large 
in comparison with the case of Vq > Eq, suggesting that the 
amplitude of oscillation does not change so rapidly. Note that 
/ is constant in the region of T < T2 both for Vq > Eq and 

For Vq > Eq, the atom feels the wide potential near the sad- 
dle point and it can oscillate with large amplitude at Ti ^ Vq. 
Since Z exhibits the peak around at T=Ti, we conclude that 
the enhancement of Z is related with the rattling effect, which 
is caused by the rapid change of the amplitude of oscilla- 
tion. For Vq < Eq, on the other hand, since the zero-point 
energy of the atom is larger than Vq, the rattling effect is 
masked by quantum fluctuations in the low-temperature re- 
gion of T < £^0- Thus, we arrive at the saturated value of Z 
dominated by zero-point oscillation at T ~ T2. In this case, 
Z does not exhibit a peak in the temperature dependence. 

Next we consider how Z is changed by 7 at T=Ti and T2 . 
For the purpose, here we define Zi and Z2 as the values of 
eq. (6) at T=Ti and T2, respectively. In Fig. 3(a), we depict 
Zi and Z2 vs. 7. For large value of 7, we find only the shoul- 
der stricture characterized by T2, as observed in the curve for 
7=0.007 in Fig. 2(a). The value of Z2 is increased with the de- 
crease of 7 and eventually, it turns to be decreased to form a 
peak at a certain value of 7, which is defined as 7*. At 7=7*, 
we also find Ti in addition to T2, as observed in the curve 
for 7=0.0005 in Fig. 2(a). When we further decrease 7 less 
than 7*, both Zi and Z2 are decreased. At very small 7, Zi 
becomes smaller than Z2, suggesting the complicated multi- 
peak structure in the temperature dependence in Z, but we do 
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Fig. 3. (Color online) (a) Zi and Z2 vs. 7 and (b) Eq and Vq vs. 7 for 



not show such a result, since it is out of the scope of this paper. 

Here we focus on the peak of Z at 7=7*. In order to clar- 
ify the meaning of 7*, we plot £^0 and Vq as functions of 7 in 
Fig. 3(b). Except for small deviation, we find that 7* indicates 
the value of 7 at which £^0 and Vq cross with each other. As 
we have emphasized above, for the case of Vb > £^0, we find 
the peak in Z atT ^ Vq due to the rattling effect. For the case 
of Vq < Eq, the ratting effect is masked by quantum fluctu- 
ations and we do not find the peak induced by rattling effect. 
At Vq ^ Eq, it is expected that both quantum and rattling ef- 
fects work cooperatively for the enhancement of Z. Thus, Z 
is considered to take the maximum for Vq Eq. 

In summary, we have discussed the electron mass enhance- 
ment factor Z in anharmonic Holstein model. It has been 
found that Z is enhanced by the rattling effect, i.e., the rapid 
change of the amplitude of oscillation. Furthermore, Z has 
been found to take the maximum value when £^0 is compara- 
ble with Vq. Namely, the heavy-electron state due to rattling 
is expected to occur significantly when zero-point energy is 
comparable with the energy at which the amplitude of oscil- 
lation is rapidly increased. 
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